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ABSTRACT 

We present an analysis of 13 of the best quality Ultraluminous X-ray source (ULX) 
datasets available from XMM-Newton European Photon Imaging Camera (EPIC) observa- 
tions. We utilise the high signal-to-noise in these ULX spectra to investigate the best descrip- 
tions of their spectral shape in the 0.3-10 keV range. Simple models of an absorbed power-law 
or multicolour disc blackbody prove inadequate at describing the spectra. Better fits are found 
using a combination of these two components, with both variants of this model - a cool (~ 0.2 
keV) disc blackbody plus hard power-law continuum, and a soft power-law continuum, dom- 
inant at low energies, plus a warm (~ 1.7 keV) disc blackbody - providing good fits to 8/13 
ULX spectra. However, by examining the data above 2 keV, we find evidence for curvature 
in the majority of datasets (8/13 with at least marginal detections), inconsistent with the dom- 
inance of a power-law in this regime. In fact, the most successful empirical description of 
the spectra proved to be a combination of a cool (~ 0.2 keV) classic blackbody spectrum, 
plus a warm disc blackbody, that fits acceptably to 10/13 ULXs. The best overall fits are pro- 
vided by a physically self-consistent accretion disc plus Comptonised corona model (DISKPN 
+ EQPAIR), which fits acceptably to 11/13 ULXs. This model provides a physical explana- 
tion for the spectral curvature, namely that it originates in an optically-thick corona, though 
the accretion disc photons seeding this corona still originate in an apparently cool disc. We 
note similarities between this fit and mode ls of Galactic bla ck hole binaries at high accre- 
tion rates, most notably the model of Done & Kubota (2005). In this scenario the inner-disc 
and corona become energetically-coupled at high accretion rates, resulting in a cooled accre- 
tion disc and optically-thick corona. We conclude that this analysis of the best spectral data 
for ULXs shows it to be plausible that the majority of the population are high accretion rate 
stellar-mass (perhaps up to 80-Af Q ) black holes, though we cannot categorically rule out the 
presence of larger, ~ 1000-M Q intermediate-mass black holes (IMBHs) in individual sources 
with the current X-ray data. 

Key words: accretion, accretion discs - black hole physics - X-rays: binaries - X-rays: 
galaxies 



1 INTRODUCTION 

Einstein X-ray observations were the first to reveal point-like, ex- 
tranuclear sources in so me nearby gala xies with luminosities in 
excess of 10 39 ergs" 1 (Fabbiano fl989h . Subsequently, many of 
these so-called Ultraluminous X-ray sources (ULXs) have dis- 
played short and long term variabili ty, which suggests they are 
predominantly accreting objects (see Miller & Colbert 2004 and 
references therein). However, the observed luminosities of most 
ULXs exceed the Eddington limit for spherical accretion onto a 
stellar-mass (~10-M©) black hole (BH). In fact, their luminosites 
are intermediate between those of normal stellar mass BH X- 
ray Binaries (BHBs) and Active Galactic Nuclei (AGN). There- 



fore, the accretion of matter onto intermediate-mass black holes 
(IMBHs, of ~10 2 -10 4 Mq) provide an attractive explanation for 
ULXs, and could represent the long sought-after 'missing link' 
between stellar mass BHs and the supermassive BHs in the nu- 
clei of galaxies. However, the large populations of ULXs asso- 
ciated with sites of act ive star formation (e.g., in the Cartwheel 
galaxy, iGaoet all [2003) demand rather too high formation rates 
of IMBHs if they are to explain the ULX class as a whole iKind 
l2004l) . An alternative to accreting IMBHs is that ULXs may b e a 
type of stellar-mass BHB with geomet rically I Ki ng et all200ll) or 
relativistically l Kordina, Falcke, & Markoff 2002) beamed emis- 
sion, such that their intrinsic X-ray luminosity does not exceed the 
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Eddington limit. Another possibility is that they are stellar-mass 
BHBs that can achiev e truly super Eddington luminositi es via slim 
lEbisa wa et aljEo03l) or radiation pressure dominated fBegelmanl 
2002) accretion discs. 

As ULXs are probably the brightest class of X-ray binary 
fueled by the accretion of matter onto a BH 1 , a knowledge of 
the properties of Galactic BHBs could be vital in interpreting 
their characteristics. Traditionally the X-ray spectra of BHBs have 
been fitted empirically with two components, namely a power- 
law continuum and a multicolour disc blackbody (MCD) compo- 
nent iMitsuda et alll984t |Makishima et al. 1986). In the standard 
picture, the power-law component is thought to represent inverse- 
Compton scattering of thermal photons from the accretion disc by 
hot electrons in a surrounding corona. As such, the power-law rep- 
resents the hard tail of the X-ray emission while the MCD compo- 
nent models the soft X-ray emission from the accretion disc. The 
MCD model itself has been formulated based on the best known 
mode l for accretion onto BHs (i.e., the thin accretion disc model, 
IShakura & Sunvae\ll973l) . 

It has long been recognised that Galactic BHBs demon- 
strate various X-ray spectral states which are defined by the bal- 
ance of these two components (i.e., power-law and MCD) at any 
one time. The three most familiar X-ray bright states are the 
low/hard (LH), high/soft (HS, also described as 'thermal domi- 
nated') and the very high (VH , or 'steep power-law') states (see 
iMcClintock & Remillardll2003l for further details). At lower mass 
accretion rates, a BHB usually enters the LH state where their X-ray 
emission is dominated by a hard power-law component (r~1.7), 
thought to arise from Comptonisation of soft photons by a hot 
optica lly thin corona 2 . I n this state, the disc is either undetected 
(e.g.. lBelloni et alll9991) or appears truncated at a much larger in- 
ner radius and hence cooler than the parameters derived for the soft 
state IWilms et alJll999llMcClintock et alj|200lh . The soft X-ray 
state is generally seen at a higher luminosity (i.e., the HS state) 
and is best explained as ~1 keV thermal emission from a multi- 
temperature accretion disc (i.e., modelled with a MCD component). 
In this state, the spectrum may also display a hard tail that con- 
tributes a small percentage of the total flux. The VH state is in many 
cases the most luminous state and is characterised by an unbroken 
power-law spectrum extending out to a few hundred keV or more. 
The photon index is typically steeper 2.5) than found in the 
LH state and generally coincides with the onset of strong X-ray 
quasi-periodic oscillations (QPOs). A MCD component may also 
be present in the VH state and the EXOSAT era demonstrated that 
some of the QPOs occur when both d isc and power-law compo- 
nents contribute substantial luminosity jvan der Klisll995l) . 

The idea of ULXs as analogues to Galactic BHBs in the 
HS state was supported by ASCA observations, which revealed 
that their 0.5-10 keV spectra were successfully fitted with the 
MCD mod el with rela tively high disc temperatures (1.0-1.8 keV, 
Makis himaetalJl2000l) . As such, the ULXs were considered to 
be mass-accreting BHs with the X-ray emission originating in an 
optically-thick accretion disc. In fact, the use of the MCD model to 
describe these spectra permits one to obtain an 'X-ray-estimated' 



1 Although the most likely reservoir of fuel for an ULX is a companion 
star, others have been suggested, for example the direct accretion of matter 
from molecular clouds i Krolik 2004). 

2 However this is still a topic of debate, with the main alternative for 
the X-ray power-law emission being synchrotr on emission from the ra- 
dio jet that is associated with t his state (e.g., iFalcke & Biermanr]|l995l 
Mark off. Falcke. & Fendej200ll) . 



BH m ass, Mxr, from the following equation (cf. Makishima et al. 
2000 equations (5)-(8)). 



(1) 



M XR = 




8.86a 



Where D is the distance to the X-ray source, which has an 
inclination i, a full bolometric luminosity (from the MCD model) 
of /boi and an observed maximum disc colour temperature 2j n . In 
addition a is the Stefan-Boltzmann constant, k is the ratio of the 
colour temperature to the effective temperature ('spectral harden- 
ing factor'), and £ is a correction factor reflecting the fact that Ti n 
occurs at a radius somewhat larger than Ri n (here, we assume that 
Ri n is at the last stable Keplerian orbit). iMakishima et alj fcOOOl) 
use values of £ = 0.412 and k = 1.7, though other work has found 
different values for the spectral hardening factor (e.g. k — 2.6 for 
GRO J1655-40. lichrader & Titarchuldl2003l) . Finally, a is a posi- 
tive parameter with a — 1 corresponding to a Schwarzschild BH. 
However, the masses inferred from the ASCA data and Equation 
(1) are far too low to be compatible with the large BH masses 
suggested by their luminosities (assuming Eddington-limited ac- 
cretio n), for stan dard accretion discs around Schwarzchild BHs. 
Makishima et al. (2000) suggested that this incompatibility could 
be explained if the BHs were in the Kerr metric (i.e., rapidly ro- 
tating objects), allowing smaller inner disc radii and hence higher 
disc temperatures. 

Chandra observ a tions have provided some support for 
the IMakishima et alj yOOOJ) results, with some ULX spec- 
tra b eing consistent with the MCD model (e.g., iRoberts et alj 
2002). However, Chandra also revealed that some ULX spec- 
tra showed a preferen ce for a power-law c ontinuum rathe r than 
the MCD model (e.e.. Istrickland et al]l200lt IRoberts et af]l2004 
iTerashima & Wilsonl l2004). It has been suggested that this pref- 
erence for a power-law spectrum could be interpreted in terms 
of the LH state seen in Galactic BHB candidates, relativisti- 
cally beamed jets or emission from a Comptonised accretion 
disc in the VH state. As well as these single component mod- 
els, ASCA and Chandra spectroscopy have also suggested the 
presence of two component spectra for some ULXs, compris- 
ing a MCD with a power-law component. For example, previous 
ASCA analyses hinted at evidence for IMBHs, i.e., cool accre- 
tion disc components (see below), but these observations were not 
sensitive enough to statistically require two component modelling 
(e.g., IColbert & Mushotzkvlll999T) . Similar results have been ob- 
tained by Chandr a e.g., ULXs in NGC 540 8 iKaaret et alj|2003l) 
and NGC 6946 jRoberts & Colbert! 120031). Conversely, Chan- 
dra spe ctra of the Antennae ULXs IZezas et all2002a : IZezas et alj 
2002b) revealed an accretion disc (MCD) component consistent 
with the high temperature ASCA results (i.e., fcTi n ~l keV), to- 
gether with a hard power-law component (r~1.2). 

It is only recently, using high quality XMM-Newton/EPIC 
spectroscopy of ULXs, that it has been demonstrated that the 
addition of a soft thermal disc component to a power-law con- 
tinuum spectrum provides a strong statistic al improvement to 
the best fitting models to UL X data (e.g., iMiller et all l2003t 
iMiller. Fabian, & Millerj l2004al) . These particular observations 
have provided strong support for the IMBH hypothesis by reveal- 
ing disc temperatures in these sources up to 10 times lower than 
commonly measured in stellar mass BHBs, consistent with the ex- 
pectation for the accretion disc around a ~1000-A/q IMBH (cf. 



It is common for the generic range of masses for IMBHs to be quoted 
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Equation (1)). However, in a few cases, XMM-Newton/EPIC ob- 
servations have revealed a more unusual two-component X-ray 
spectr um. A detailed analysis of su ch a source is presented in 
Stobb art. Roberts. & Warwick! fe004l) . In this case, the ULX is the 
brightest X-ray source in the nearby (1.78 Mpc) Magellanic-type 
galaxy NGC 55 and, although its X-ray luminosity only marginally 
exceeds 10 39 erg s^ 1 , it represents one of the highest quality ULX 
datasets obtained to date. The initial low flux state data were best 
fitted with an absorbed power-law continuum (r~4), while a sub- 
sequent flux increase was almost entirely due to an additional con- 
tribution at energies > 1 keV, adequately modelled by a MCD com- 
ponent (fcT in ~0.9 keV). 

Whilst this accretion disc component is reasonable for stellar- 
mass BHs, the dominance of the power-law continuum at soft X- 
ray energies is problematic. Such a soft power-law cannot repre- 
sent Comptonised emission from a hot corona, as one would not 
expect to see the coronal component extend down below the peak 
emissivity of the accretion disc, where there would be insufficient 
photons to seed the corona. Alternative sources of seed photons for 
the corona are unlikely; for example, the incident photon flux of the 
secondary star at the inner region s of the acc retion disc is too low to 
provide the seeding (cf. lRoberts et all2005l) . It also seems unlikely 
that the power-law emission could arise from processes at the base 
of a jet, as these are typically represented by much harder photon 
indices than measured here (F~ 1.5-2; Ma rkoff et all20 05 and ref- 
eren ces therein). I ndeed, with the possible exception of NGC 5408 
X-l iKaaret et all2003l) . there is no evidence that ULXs do display 
bright r adio jets, though this cannot be excluded by current obser- 
vations l Kordina et al. 2005). The possibility of the soft component 
resulting from an outflow of material from the accretion disc may 
also be discounted as this would produce a thermal spectrum rather 
than a power-law continuum. 

Although this spectral description has not been seen in Galac- 
tic systems, a second case has been reported independently for the 
neares t persistent extragalactic ULX (M33 X-8) by Foschini et al. 
(2004). The non-standard model provided the best fit to this ULX 
with r~2.5 and fcTi n ~1.2 keV. However, this source is also at the 
low luminosity end of the ULX regime with Lx ~2 x 10 39 erg s~ 1 . 
A further possible case, in a more luminous ULX, h as arisen from 
the X MM-Newton data analysis of NGC 5204 X-l jRoberts etail 
120051) . In this case the authors show that there is spectral ambi- 
guity between the non-standard fit (r~3.3, fcTi n ~ 2.2 keV) and 
the EVIBH model (T~2.0, fcT in ~ 0.2 keV), with both providing 
statistically acceptable fits to the data. Even more recently, two ad- 
ditional examples of this spectral form have been unco vered in an 
XMM-Newton survey of ULXs bv lFeng & Kaaretl <2005h . 

Although it is difficult to derive a literal physical interpretation 
from the non-standard model, it does provide an accurate empirical 
description of ULX spectra in some cases, and as such it has the 
potential to provide new insights into the nature of these sources. 
Therefore in this paper we re-evaluate current data in an attempt to 
determine the best spectral description for the shape of high qual- 
ity ULX spectra, and ask what consequences this has for the idea 

as 20-10 6 Mq . The lower limit comes from a c onsideratio n of the mea - 
sured masses of BHs in our own Galaxy <McClintock & RemiUarJl2003l) . 
and a theoretical limit for the mass of a BH formed from a single massive 
star iFrver & Kaloaera 2001). However, more recent population synthesis 
analyses show that BHs of up to ~ 80- Mp, may be formed in young stel- 
lar populations 1 Belczvnski et al. 2004). Hence, when we refer to IMBHs 
in this paper we refer specifically to the larger ~1000-Mq IMBHs implied 
by the cool accretion disc measurements. 



of ULXs as accreting IMBHs. The paper is structured as follows: 
Sec. [5]- introduction to the ULX sample; Sec. [3]- details of the 
observations and data reduction; Sec. [4]- description of the spec- 
tral analysis; Sec. [5]- a comment on the luminosity and inner disc 
temperature relationship of these ULXs; Sec. [6]- a discussion of 
our results; and finally Sec.Q- our conclusions. 



2 THE SAMPLE 

As our primary goal is to find the best description(s) of the shape 
of ULX spectra, only the highest quality datasets were chosen. 
The ULXs were initially s elec ted from the ROSAT c atalogues of 
iRoberts & Warwickl feOOfJ) and lColbert & Ptakl (2002) to provide a 
list of historically-bright ULXs that are resolved at a spatial resolu- 
tion similar to XMM-Newton 4 . We applied a source selection crite- 
ria of observed count rates of > 10 counts ks -1 in the ROSAT KRI 
camera, combined with > 10 ks of XMM-Newton/EPIC data avail- 
able in the archive by December 2004, to select ten ULXs with 
potentially sufficient counts for very detailed spectral analysis. In 
addition, we included three more high quality ULX datasets: pro- 
prietary data for Holmberg II X-l (hereafter Ho II X-l), and two 
sources not quite bright enough in the ROSAT bandpass to be clas- 
sified as ULXs, namely the ULX in NGC 55 and M33 X-8. Whilst 
some of the sources in this sample have been observed more than 
once, we have only selected the longest individual exposure in each 
case to provide the clearest single view of their spectra. Our fi- 
nal sample of 13 sources from 12 different galaxies is listed in Ta- 
ble □ The selected ULXs are located at distances of between 800 
kpc and 17.8 Mpc, possess XMM-Newton count rates between 0.1 
and 8.9 count s _1 and cover the full range of ULX luminosities 
(~10 39 erg s _1 - few xlO 40 erg s _1 ). Hereafter we refer to the 
sources by their names as given in column (1) of TableQ 

3 OBSERVATIONS AND DATA ANALYSIS 

In this work we hav e utilised data from the EP IC cameras on 
board XMM-Newton <Turner et alJl200ltlStruder et alj|200ll) . The 
datasets were obtained through the XMM-Newton public data 
archive (excluding proprietary Ho II X-l data) and details of the 
observations are shown in Table [2] The data were processed and 
reduced using the standard tools of XMM-SAS software v.6.0.0. In 
some cases the observations were affected by soft proton flaring for 
which preliminary cleaning was necessary. For these observations 
we extracted full field X-ray (0.3-10 keV) light curves and screened 
for flaring using Good Time Interval (GTI) files based on either a 
time or count rate criterion. The NGC 55 ULX and NGC 5204 X-l 
observations were not affected by flaring episodes. The NGC 2403 
X-l observation was only affected by flaring at the end of the expo- 
sure (the last ~20 ks), therefore we used a time selection to exclude 
this flaring event. For the remaining sources we used a count rate 
cut-off criterion to produce a GTI file. The exact value of the cut-off 
was allowed to vary from field-to-field, to provide the best compro- 
mise in each case between excluding high background periods and 
facilitating the longest available exposure on the ULX. In practise 
the actual cut-off values varied in the 6.5-17.5 count s _1 (0.3-10 
keV) range. In all cases we used only those data recorded when the 

4 The ROSAT lists were deemed appropriate as most bright ULXs are per- 
sistent a nd vary by factors of no more than 2-3 in flux over a baseline of 
years, cf. Ro berts et af]<2004 . 
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Table 1. The sample 











7V H 


d 


Lx 


Source 


Alternate names 


K.A. (J2UUU) 


Ub,L. (JzUUU) 


(1U cm ) 


(Mpc) 


(1U erg 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


NGC 55 ULX 1 


XMMU J001528.9— 391319 a 
NGC 55 6 c 
Source 7 d 


00 15 28.9 


-39 13 19.1 a 


1.74 


1.78 6 


1.3 


M33 X-8 2 




01 33 50.9 


+30 39 37.2 e 


5.69 


0.70 f 


1.0 


NGC 1313 X-l 3 


Source 6 9 
IXO 7 i 


03 18 20.0 


-66 29 11.0 h 


3.96 


3.70" 


4.7 


NGC 1313 X-2 3 


Source 4 9 
IXO 8 i 


03 18 22.3 


-66 36 03.8 k 


3.94 


3.70' 


1.7 


NGC 2403 X-l 4 


Source 21 1 


07 36 25.5 


+65 35 40.0 1 


4.17 


4.20 


2.7 


Ho II X-l 5 


1X0 31 i 


08 19 29.0 


+70 42 19.3 m 


3.41 


4.50 1 


17 


M81X-9 6 


HolmbergIXX-1 n 
IXO 34 i 
NGC 3031 10 i 
H44 r 


09 57 53.2 


+69 03 48.3 ° 


4.25 


3.55 p 


12 


NGC 3628 X-l 4 


IXO 39 J 


11 20 15.8 


+ 13 35 13.6 s 


2.22 


7.70 / 


5.2 


NGC 4395 X-l 4 


NGC 4395 X2 4 
IXO 53 3 


12 26 01.5 


+33 31 30.5 " 


1.36 


3.60 / 


0.6 


NGC 4559 X-l 4 


X-7 " 
IXO 65 i 


12 35 51.7 


+27 56 04.1 w 


0.82 


9.70 ^ 


9.1 


NGC 4861 ULX 1 


IXO 73 j 
XI 9 


12 59 01.9 


+34 51 13.5 x 


1.21 


17.80 i 


8.8 


NGC 5204 X-l 4 


IXO 77 i 
HST3 z 

Ul aa 


13 29 38.6 


+58 25 05.7 z 


1.38 


4.80 f 


4.4 


M83 ULX 1 


Source 13 bb 
IXO 82 » 
H30 dd 
H2 ee 


13 37 19.8 


-29 53 48.9 cc 


3.69 


4.70 1 


1.0 



NOTES: (1) Source designation; (2) Alternative names; (3^1) X-ray source position from XMM-Newton and Chandra data, or posi- 
tion of p ossible optical coun terpart; (5) Galactic absorption column density from the 'NH' FTOOLS program (based on the measure- 
ments of iDickev & Loc kman 1990); (6) Distance to the host galaxy; (7) Observed X- ray luminosity (0.3-10 k eV ) based on the result s 
o f the physically self- c onsist ent modelling (see l a ter). R EF ERENCES: 1 This pa per. 2 | M arkert & RalliJ Il983l), 3 IC olbert et all Il995t) . 
4 iRoberts & Warwick! f200(l). 5 I5ewangan et al] 120041). 6 iFabbianol Jl983). " IStobbart et alj l2004l). b iKarachentsev et" 2003 1, 
[ Read. Ponman & Strickland! Il997t). d [Schlegel. Barrett & Singhl j 199711. 
Isch] 



Read. Ponman & Stricklan 

h 



nd 
M 



1997F7 " ISchleeel. Barrett & Sing! 
Miller et alj l2003l> . i iTullvl ll98j 



irkert & R 
etalJ 120, 

19971). e iFoschini et all |2004I).^Tho. Filippenko. & Sargenll I 

JColbert & Ptakl 12002I). k IZampieri et alj |2004 ' ISchleeel & Pamiutil 

iMiller et all l2004al). iRamsev et alj <2006l). p IPaturel et all l2002t) . 9 iRadeckei ll997t) . 



1997) 



Schleg el et alj 120001) 
120031). m iKaaret Ward. & Zeza 

r llmmler & Wangj l200ll) . s IStrickland et alj 12001ft. f iLira. Lawrence. & Johnso: 
dra X-ray Centre: http://cxc.harvard.edu/chaser " Vogler. Pietsch. & Bertoldi 1 199 

0141 15010) via XMM-Newton Science Archive: http://xmm.vilspa.esa.es/external/xmm_data_acc/xsa/index.shtml y Liu & Bregman 1 20 



i20()0), " Source-l i st (Ob s.Id 5302) via Chan- 
. w ICropper et alj <2004. x Sourcelist (Obs.Id 



iGoad et ail <200_i). aa Liu. Bregman. & SeitzeJ <2004l 
iTrinchieri. Fabbiano. & Paulumbd<1985l) . 



Ehleetal. 1 1998), cc Soria & Wu (2002), dd Immleretal. 



MOS and pn cameras were in simultaneous operation (see Table|2| 
for the net good exposure times). 

For each ULX, events were extracted in a circular aperture 
centred on the source position given in Table Q The background 
was taken from a circular region, near to the source in the pn camera 
and at the same distance from the readout node. The chosen back- 
ground regions were the same in all three detectors where obser- 
vations were taken using the full-frame observing mode. However, 
in three cases the MOS cameras were operated in either the large- 
window or small-window mode, so here we used a background re- 
gion closest to the one used for the pn extraction. The size of the 
source and background spectral extraction regions are listed in Ta- 
ble^ 

For the pn camera the spectra were extracted using event pat- 
terns 0-4, which allows 'single' and 'double' pixel events. We also 
set 'FLAG=0' to exclude all events at the edge of the CCD and 



events from bad pixels. For the MOS detectors we used event pat- 
terns 0-12 which allows 'single', 'double', 'triple' and 'quadruple' 
pixel events. We employed a less conservative screening criterion 
for the MOS data by using the flag expression #XMMEA_EM to ex- 
clude hot pixels and events outside of the field of view. The SAS 
task especget was used to produce source and background spectra 
for each ULX, together with the appropriate Redistribution Ma- 
trix File (RMF) and Ancilliary Response File (ARF). Spectral files 
were grouped to require at least 20 counts bin -1 before fitting, to 
ensure Gaussian statistics. 

Model spectra were fitted to the data using the HEAsoft X- 
ray spectral fitting package XSPEC (v. 11. 3.0). The pn, MOS-1 and 
MOS-2 spectra for each object were fitted simultaneously, but we 
included constant multiplicative factors in each model to allow for 
calibration differences between the cameras. This value was frozen 
at unity for the pn data and allowed to vary for the MOS detectors, 
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Table 2. XMM-Newton Observation Log 



Source Obs.ID Date Duration (s) Net exp. (s) Rate(cts _1 ) MOS-1 MOS-2 pn Position Ref 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



NGC 55 ULX 


0028740201 


2001-11-14 


34025 


30410 


2.08 


FF 


FF 


FF 


on-axis 


1 


M33 X-8 


0102640101 


2000-08-04 


18672 


6850 


8.55 


sw 


SW 


FF 


on-axis 


2,3,4 


NGC 1313 X-l 


0106860101 


2000-10-17 


42769 


18490 


1.15 


FF 


FF 


FF 


on-axis 


4,5,6 


NGC 1313 X-2 


0106860101 


2000-10-17 


42769 


18490 


0.39 


FF 


FF 


FF 


off-axis™ 


4,5,6 


NGC 2403 X-l 


0164560901 


2004-09-12 


84600 


57063 


0.47 


FF 


FF 


FF 


on-axis 


4 


Ho II X-l 


0200470101 


2004-04-15 


111999 


47260 


4.52 


LW 


LW 


FF 


on-axis 


7 


M81 X-9 


0112521101 


2002-04-16 


11935 


8440 


3.25 


FF 


FF 


FF 


on-axis 


4.5 


NGC 3628 X-l 


0110980101 


2000-11-27 


60745 


45260 


0.20 


FF 


FF 


FF 


on-axis 


4.5 


NGC 4395 X-l 


0142830101 


2003-11-30 


118900 


98900 


0.18 


FF 


FF 


FF 


on-axis 


4.8 


NGC 4559 X-l 


0152170501 


2003-05-27 


43290 


37660 


0.48 


SW 


SW 


FF 


on-axis 


4.9 


NGC 4861 ULX 


0141150101 


2003-06-14 


29600 


14590 


0.11 


FF 


FF 


FF 


on-axis 




NGC 5204 X-l 


0142770101 


2003-01-06 


32999 


17048 


0.87 


FF 


FF 


FF 


on-axis 


4,10 


M83 ULX 


0110910201 


2003-01-27 


31722 


21130 


0.21 


FF 


FF 


EFF 


off-axis' 


4 



NOTES: (1) Source designation; (2) Observation identifier; (3) Observation date (yyyy-mm-dd); (4) Observation duration (NB. this does include cali- 
bration observations in some cases); (5) Useful exposure after correcting for flaring episodes and ensuring simultaneous operation of the EPIC cameras; 
(6) Combined EPIC count rates derived from the combined X-ray light curves (0.3-10 keV); (7-9) Observing mode of each EPIC detector (SW : Small 
Window, LW : Large Window, FF : Full Frame, EFF : Extended Full Frame); (10) Source position with respect to the centre of the pn field of view: 
a— ~7 ' o ffset; b- ~6.5 ' offse t; (1 1) References for previou s analyses of these dat asets (although not all referenc es contain a c om plete analysis of th e 
U LX): 1-IStobbart et alJfcoOft. 2-lFoschini et alJ<2004l). 3jPietsch et al] <2004l). 4-lFeng & KaareJl2005t) . 5- IWang et aill2004l) . 6 4Miller et aT]l2003l) . 
7 4Goad et alj<2006l) . 8- IVaughan et alj|2005t) . 9 JCropper et alj<2004l) . 10 4Roberts et alj|2005t) . 



with the values typically agreeing within 20 per cent. Quoted fluxes 
are an average of the three measurements. 

Spectra were initially fitted in the 0.3-10 keV band. However, 
in several cases we found that the data below 0.5 keV did not agree 
with any of the models we tried, leaving large residuals (particu- 
larly with respect to the pn data) to the best fits. We note that the 
datasets with this problem - NGC 55 ULX, NGC 1313 X-l & X-2 
and NGC 3628 X-l - were all obtained in the years 2000-2001 (in 
fact M33 X-8 is the only observation from this epoch without these 
fitting residuals). At first glance this would appear to be a problem 
with the calibration at early mission times; however the pn cali- 
bration at least has remained remarkably static in orbit (EPIC team, 
priv. comm.). As we cannot explain this variation we conservatively 
excluded these data in all EPIC cameras, such that spectral fitting 
was restricted to the 0.5-10 keV range in these four cases. Also, in 
the case of Ho II X-l there is a possible pn calibration feature that 
is particularly prominent at low energies I Goad et al. 2006), there- 
fore we followed the method of these authors and excluded the pn 
data below 0.7 keV to account for this, whilst retaining MOS data 
down to 0.3 keV. 

The X-ray spectra were modified for absorption follow- 
ing B alucinska-Church & McCammonl | l992l) . assuming the solar 
abundances of lAnde^^^Grevesseni989l) (i.e.. the WABS compo- 
nent in XSPEC). We used two absorption components, one of which 
was fixed for each source to represent the appropriate foregroun d 
column density through our Galaxy I Dick ey & Lockmar] [l990 - 
see Table0, and the second component was left free to fit the data 
to represent additional absorption within the host galaxy and/or in- 
trinsic to the ULX. The errors quoted in this work are at the 90% 
confidence level for one interesting parameter. Throughout this 
analysis we distinguish statistically acceptable fits from unaccept- 
able fits using a fixed criterion of P rc j < 95 per cent, where P ro j 
is the probability of rejection derived directly from the x' 2 statistic 
for the spectral fit. 



Table 3. Spectral extraction apertures 



Extraction radius 
Source Source Background 



(1) 


(2) 


(3) 


NGC 55 ULX 


60 


75 


M33 X-8 


35 


120 


NGC 1313 X-l 


35 


70 


NGC 1313 X-2 


35 


70 


NGC 2403 X-l 


45 


45 


Ho II X-l 


45 


60 


M81 X-9 


30 


45 


NGC 3628 X-l 


35 


70 


NGC 4395 X-l 


40 


40 


NGC 4559 X-l 


40 


80 


NGC 4861 ULX 


15 


45 


NGC 5204 X-l 


40 


50 


M83 ULX 


35 


50 


NOTES: (1) Source 


designation; 


(2)-(3) 



Aperture radii in arcseconds 

4 SPECTRAL PROPERTIES 
4.1 Simple models 

The X-ray spectra appeared relatively smooth and featureless in 
each case, hence we began by fitting simple continuum models to 
the data. We tried the two models most often used in the past as 
single-component descriptions of ULX spectra, namely a power- 
law continuum model and the canonical MCD model. The details of 
these fits are given in Tableland the quality of the X-ray spectra is 
illustrated using the power-law fits in FigureQ Due to the excellent 
spectral quality, we were able to reject the power-law model in 8 
out of the 13 cases at > 95% confidence. Interestingly, the five 
ULXs for which a power-law is an adequate description of the data 
(NGC 1313 X-l & X-2, M83 ULX, NGC 4861 ULX, NGC 5204 
X-l) include the three poorest quality datasets. Even more notably, 
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Figure 1 - continued 



the MCD model does not provide an acceptable fit to any of the 
spectra. 



4.2 Power-law + MCD model 

As simple spectral models were inadequate, we proceeded to fit the 
data with the two component model generally employed for BHB 
systems, namely the combination of a power-law plus a MCD com- 
ponent. More specifically, we began by testing the standard IMBH 
model, i.e., a cool accretion disc plus power-law continuum model. 
This spectral description improved on the simple models by provid- 
ing acceptable fits to 8 sources (NGC 1313 X-l & X-2, M81 X-9, 
NGC 3628 X-l, NGC 4559 X-l, NGC 4861 ULX, NGC 5204 X-l, 
M83 ULX), with 0.1 keV < kT in < 0.3 keV and 1.6 < T < 2.5 



(Table[5}- As with previously published work, these disc tempera- 
tures are broadly consistent with IMBHs of around -^lOOO-M© in 
size, though the power-law slopes are puzzlingly shall ow for w hat 
are supposed HS (or VH state) sources (see Robe rts et alJ 2005 for 
further discussion). In all cases the 0.3-10 keV flux is dominated 
by the power-law component. 

We then attempted to fit the X-ray spectra with the non- 
standard model i.e., with the power-law component dominant at 
soft energies. This also provided statistically-acceptable fits to 8 
ULXs (M33 X-8, NGC 1313 X-2, NGC 2403 X-l, NGC 3628 X- 
1, NGC 4559 X-l, NGC 4861 ULX, NGC 5204 X-l, M83 ULX), 
with disc temperatures of 1.2 keV < kT[ n < 2.2 keV (excluding 
M 83 for which the temperature is unconstrained below ~4 keV) 
and 2.5 < F < 4.5 (except NGC 3628 X-l, for which F is un- 
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Table 4. Single component spectral fits 



WA WA rU 




1 


j\p 


X /ao1 


NGC 55 ULX 


3.50±0.08 


3.38±0.03 


3.09+0.08 


1128.9/830 


M33 X-8 


1.94±0.06 


2.28±0.02 


6.02+0.12 


1790.1/1115 


NGC 1313 X-l 


1.05±0.11 


1 84+ ' 04 
'■^-o.os 


0.55+0.02 


723.0/675 


NGC 1313 X-2 


1.68±0.02 


9 97+0.08 
z - z '-0.07 


0.31+0.02 


261.7/257 


NGC 2403 X-l 


3.78±0.13 


2.40±0.03 


0.65+0.02 


1130.2/851 


Ho II X-l 


1 53+0 04 


2 76+0 02 


3 23+0 04 


1690 6/1314 


M81 X-9 


1.57±0.08 


1.89+0.03 


, O9+0.06 


938.1/864 


NGC 3628 X-l 


2.77±0.03 


1.45+0.05 


0.10+0.01 


459.8/406 


NGC 4395 X-l 
NGC 4559 X-l 


~ _~-i-0.17 
*• IJ -0.16 

1.15±0.08 


- T-,+0.12 

4 - 37 -o.n 

o to + 0.05 
2 - JK -0.04 


0.09+0.01 
0.27+0.01 


540.8/368 
752/599 


NGC 4861 ULX 
NGC 5204 X-l 
M83 ULX 


. « A -4-0 31 

0.50±0.08 
i 09+0.21 
^-0.20 


z -^ y -0.15 
2 iQ+0.05 
zlu -0.04 
2 ci+0.11 

Z - J -0.10 


0.09+0.01 
19+0.02 

0.14+0.01 


94.0/74 
572.3/529 
199.2/207 


W A * WA * D I S KB B 




fcii n 


a e 

AmCD 


^,2 lAnf 

X /dot 


NGC 55 ULX 


0.37 


0.60 


1.62 


1796.3/830 


M33 X-8 





1.08+0.01 


62+ 02 

u - oz -0.01 


1538.7/1115 


NGC 1313 X-l 





1.32 


0.04 


1409.6/675 


NGC 1313 X-2 


< 0.02 


0.93+0.04 


0.05+0.01 


443.3/257 


NGC 2403 X-l 


, 20+0.08 


1.04+0.02 


0.06+0.01 


923.7/851 


Hon x-i 





0.59 


2.74 


8646.7/131 


A/IQ 1 V Q 


n 
U 


1.42 


0.09 


1 / JO.0/O04 


NGC 3628 X-l 
NGC 4395 X-l 


1.04±0.02 
0.13±0.09 


9 9 7 +0.13 
z - z '-0.12 

0.32+0.01 


(1.38+°-^) x 10 
0.8li°;« 


" 3 466.6/406 
693.0/368 


NGC 4559 X-l 





0.69 


0.15 


1651.4/599 


NGC 4861 ULX 





0.64 


0.06 


166.9/74 


NGC 5204 X-l 





0.66 


0.33 


1650.1/529 


M83 ULX 


< 0.01 


0.64+0.03 


0.11+0.02 


310.7/207 


NOTES: Models are 


abbreviated 


to XSPEC syntax: WA-absorption 


components for the 



Galactic value and external absorption; PO-power-law continuum; DISKBB— MCD. a Ex- 
ternal absorption column (10 21 atoms cm" 2 ), b PO photon index, c PO normalisation 
(10 -3 photon cm" 2 s" 1 keV -1 at 1 keV ), d Inner disc temperature (keV). e DISKBB 
normalisation (((Ri n /km)/(D/10 kpc)) 2 cos ©; where Ri n -inner disc radius, D-distance to 
source, 0-inclination angle of the disc). Where the reduced \ 2 value (i.e. ¥y) exceeded 
2, errors are not shown. Statistically acceptable fits are highlighted in bold. 



constrained below ~1.2) (Table|5J. The balance in 0.3-10 keV flux 
between the two components was much more varied for this model, 
though most sources showed a relatively equitable balance (ratios 
of no more than 3 : 1 in either direction, and ~ 1 : 1 in five cases) . 

In total, the power-law + MCD combination could not be re- 
jected at > 95 per cent confidence for 10/13 sources in the sam- 
ple. Two of these sources were unambiguously best fitted with the 
IMBH model (NGC 1313 X-l, M81 X-9) while two more were 
best fitted by the non-standard model (M33 X-8, NGC 2403 X-l). 
Statistically acceptable fits existed for both models for the remain- 
ing six sources. Three sources in this sample (NGC 55 ULX, Ho II 
X-l and NGC 4395) rejected both models at > 95 percent confi- 
dence, although a comparison of the fits show that the spectrum of 
the NGC 55 ULX appeared to fit much better to the non-standard 
model, whilst Ho II X-l and NGC 4395 X-l much preferred the 
IMBH model. 

This sample of high quality ULX datasets demonstrates spec- 
tral ambiguity for six sources (NGC 1313 X-2, NGC 3628 X-l, 
N GC 4559 X-l NGC 4861 ULX, NGC 5204 X-l, M83 ULX). 
As iRoberts et alJ j2005l) suggest, a key discriminator between the 
models could be a test for spectral curvature at high energies (> 2 
keV). Curvature is not expected in the IMBH model where the 



power-law component is dominant above 2 keV, whereas sources 
best-fitted by the non-standard model have the curved MCD com- 
ponent dominant at these hard X-ray energies. We therefore at- 
tempted to fit the high energy (2-10 keV) ULX spectra with a 
broken power-law model. For comparison we also fitted the 2-10 
keV ULX spectra with a single power-law component, as shown 
in Table [6] We did not include absorption in the fit in either case, 
as single- and two-component fits generally limit absorption to 
< 4 x 10 21 atoms cm -2 , which should not strongly affect the 
data above 2 keV. When we do include absorption in these fits 
we actually measure far higher columns in several cases, which is 
physically unrealistic (and may be a result of the absorption com- 
pensating for intrinsic curvature). The results of the broken power- 
law spectral fits, together with the statistical probability of the fit 
improvement over the single power-law fits (using the F-test) are 
shown in TableQ 

To demonstrate the validity of this method, the two sources 
with unambiguous non-standard fits (M33 X-8 and NGC 2403 X- 
1), plus NGC 55 ULX which clearly prefers this model, show unac- 
ceptable power-law fits that are made acceptable by the inclusion of 
a break. In these cases the improvement is highly statistically sig- 
nificant (> 9<t improvement over a simple power-law fit accord- 
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Table 5. Two component spectral fits 



WA*WA*(PO+DISKBB) 1 



TV- 



II 



kT in 



nicD 



: (xl0 2 ) 



X 2 /dof 



/ 



NGC 55 ULX 
M33 X-8 
NGC 1313 X-l 
NGC 1313 X-2 
NGC 2403 X-l 
Ho II X-l 
M81 X-9 
NGC 3628 X-l 
NGC 4395 X-l 
NGC 4559 X-l 
NGC 4861 ULX 
NGC 5204 X-l 
M83 ULX 



+0.42 
0.09 
+0.25 
0.06 
+0.52 
0.48 
+0.76 
0.42 
0.04 
0.05 
+0.10 
0.09 
+0.43 
0.23 
42 g+0.77 
^■ ZD -0.70 

1 96+ ' 26 

1 - yD -0.40 
7 ,0+0.70 
" J -0.20 

1 75+ 1 13 

' -0.69 

0.68+ ' 22 



3.50 
4.43 
2.37 
1.80 
7.18 
1.74 
2.17 



1.28 



0.18 
+0.54 
0.37 



3.38±0.03 
7 s-,+0.02 

z - JZ -o.oi 
1.77±0.06 

7 I ,+0.15 
zal -0.17 
- fi o+0.02 
2 - 68 -0.05 

2.63±0.03 

, O1+0.05 
'• S1 -0.04 

, t;o+0.07 

L5S -0.04 
,+0.18 
-0.27 
,+0.07 
-0.05 
2 24+0-23 
Z ' Z +-0.24 

1 qi+0.07 

2 4 7 +0.14 
z -^'-0.17 



□+0.07 
-0.08 



3.72] 

2.23" 1 



3.09] 

,+0.20 

-0.12 

0.51±0.04 

; +0.05 
-0.06 

0.98±0.03 
2.81±0.09 

i+0.14 
-0.10 

0.12±0.01 
0.05±0.01 

,+0.03 
-0.01 

0.07±0.02 
0.31±0.03 
14+ 002 

u -'^-0.03 



8.5 r 



0.25"* 



1.74 T 



0.24"* 



< 0.16 

0.09±0.002 

0.19±0.02 
027 +o.io 

u - z '-0.09 
g+0.002 
"■"^-0.004 

0.19±0.01 
0.20±0.03 
0.08±0.01 
0.18±0.02 
0.14±0.Q1 
0.18 
0.22] 
[~ 0.2] ■ 



< 9.71 x 10 4 

hl694 
-1055 
,+1.39 
-0.55 
j+0.27 



33127 

0.81" 1 



,+0.08 
-0.04 
,+0.04 
-0.03 
k 



0.031 002 

+644 
JZ '-224 
,+1.23 
-0.78 

1 05+ 2 - 78 
'■UJ-0.61 

7+402 

-96 

1+0.09 

-0.05 

rf+5.64 

-0.96 

,+0.21 
-0.07 



2ir 



127] 
O.llj 
2.74" 1 



0.12] 
< 0.98 



1128.8/828 
1527.0/1113 
664.5/673 
255.0/255 
1028.5/849 
1453.1/1312 
877.2/862 
424.6/404 
492.2/366 
597.2/597 
84.4/72 
529.9/527 
198.6/205 



99.9 
92.7 
90.9 
91.7 
93.6 
89.7 
94.0 
98.1 
71.6 
81.5 
82.4 
88.2 
96.6 



NGC 55 ULX 
M33 X-8 
NGC 1313 X-l 
NGC 1313 X-2 
NGC 2403 X-l 
Ho II X-l 
M81 X-9 
NGC 3628 X-l 
NGC 4395 X-l 
NGC 4559 X-l 
NGC 4861 ULX 
NGC 5204 X-l 
M83 ULX 



1.43" 1 



0.60^ 



H0.16 
-0.17 

16+ 010 

ulo -0.09 

1+0.14 
-0.25 

< 0.54 
2.02±0.30 

< 0.01 

64+ 14 
u - 04 -o.io 

7 ,+0.28 
u ' /1 -0.24 

0.34+™ 

0.801°^ 
4 6 +o.05 

u>to -0.03 

< 0.03 

< 0.37 



0.20±0.01 

0.27±0.02 
25 +0.02 
u - ZJ -o.oi 

u - z '-0.02 
02 ,+0.03 
u ' zl -0.02 

0.23±0.002 
0.27±0.01 
60+ 05 

u - DU -0.06 

0.17±0.01 

0.17±0.01 

0.20±0.03 

0.20±0.01 
021 +0.02 

u ' zl -0.03 



O.J 

o.4 r 



, ,5+0.33 
-'■-'-'-0.31 

3 09+ 31 
6 +o.li 

J -0.10 
,+0.07 
-0.06 
n ,o+0.14 
U - JS -0.09 

4 ,7+0.06 
^■-"-0.07 

2.56±0.16 
16+ ' 04 

u - lo -0.06 

0.15±0.02 
,+o.n 

-0.08 

;+0.09 

-0.04 

0.67±0.03 

t;+0.05 

-0.02 



0.61] 
0.15" 1 



0.16"* 



0. 81±0.01 
1 26+ 02 

1, ZD -0.03 

2.20±0.10 

,+0.15 
-0.16 
,+0.04 
-0.03 
1.28±0.02 
,+0.12 
-0.11 
,+0.74 
-0.87 

0.61±0.05 

, ,,+0.07 
'"' -0.06 
, ,7+0.20 
1 -"_0.08 

1 69+° 10 

'• Dy -0.09 
-,+0.13 
-0.05 



1.71] 
1.12"* 



2.21] 
4.32"* 



1.10"* 



H53 
-47 

303 ±?6 
7+0.31 
-0.16 
,+0.49 
-1.36 

46.63±3.42 

77.69l«; 44 

15.93+ 3 ; 34 

13+0- 17 
u- 1 -0.07 
-;+9.56 
-6.26 
,+2.29 
-1.82 
1+1.89 
-1.13 
7 +1.69 
-1.40 
,,+2.76 
-3.55 



343" 1 



5.07] 
4.02" 1 



17.66]* 
9.94+ 
2.38j 
7.37] 1 
9.28"* 



884.1/828 
1216.7/1113 
666.5/673 
257.1/255 
830.0/849 
1677.9/1312 
922.3/862 
429.4/404 
464.5/366 
604.9/597 
77.7/72 
527.8/527 
193.2/205 



36.4 
11.7 
17.7 
27.7 
8.7 
41.4 
17.5 
17.7 
65.5 
30.6 
34.1 
28.0 
27.8 



0.1 

7.3 

9.1 

8.3 

6.4 

10.3 

6.0 

I. 9 
28.4 
18.5 
17.6 

II. 8 
3.4 



WA*WA*(PO+DISKBB) 2 




T b 




kT in d 


A M CD e (xl0- 3 ) 


X 2 /dof 


fv f 


X MCD 


NGC 55 ULX 
M33 X-8 


4 32+0.24 
^■ JZ -0.34 
1 42+ ' 21 


4 31+0-19 
H - J -0.27 
240+O.I8 
z -^-0.14 


3 7Q+0.21 
-0.27 

J - 1J -0.32 


0.86±0.03 
1.18±0.06 


170 47+ 35 - 51 
264 5<;+53.06 


957.6/828 
1187.0/1113 


67.6 
48.8 


32.4 
51.2 


NGC 1313 X-l 


















NGC 1313 X-2 
NGC 2403 X-l 
Ho II X-l 


, ,,+0.97 
z - JZ -0.62 
4 77+1-01 
H -"-0.99 
2 34+0.13 
Z - J -0.12 


7 q,+0.75 
z '"-0.47 

4 05+0.70 

345+O.II 

J - H -o.io 


36+ 09 

U - JO -0.05 

55+ - 20 

"""-0.15 
, 70+O.H 

j - /s -o.io 


2 65+1-06 

z - DJ -0.78 

1 16+ ' 04 
'■'°-0.05 
1 70+O.IO 

^'^-o.oo 


<2.19 
36 93+10.14 

10.33+. 3 ; 33 , 


255.9/255 
847.7/849 

1508.8/1312 


64.4 

27.3 
74.5 


35.6 
72.7 
25.5 


M81 X-9 


















NGC 3628 X-l 


1 43+ ' 47 
1 - HJ -0.26 


< 1.18 


< 0.04 


1.40±0.21 


4.68±^ 


439.5/404 


53.4 


46.6 


NGC 4395 X-l 


















NGC 4559 X-l 
NGC 4861 ULX 
NGC 5204 X-l 


, ,c+0.46 
J - J -0.41 
, ,, + 1 .75 
-'■-'-'-1.30 

1 64+ ' 37 


4 47+0-37 
^'-0.34 

4 36+1.34 

^" ,o -1.08 
, ,,+0.34 
-'"'-'-0.32 


41+0 05 

-0.04 
13+0- 08 

44+0-05 
u -^*-0.03 


1 61+ 013 
1 60+"' 60 

'■ DU -0.32 
>y -1-0.32 
z - z -"-0.26 


4 04+1-64 
^■"^-1.59 

1 09+ 1 ' 34 

'• uy -0.87 

1 96+ 117 

1 - yD -0.87 


637.0/597 
85.6/72 
525.5/527 


47.3 
51.9 
50.3 


52.7 
48.1 
49.7 


M83ULX 


, ,7+0.92 
^■"-o^s 


7 7 t:+0.94 
z -' J -0.45 


0.14±0.02 


< 4.03 


<4.53 


198.1/205 


82.9 


17.1 


WA*WA*(BB+DISKBB) 




kT h 




kT in d 


A M CD e (xl0- 3 ) 


X 2 /dof 


IV J 


'^MCD 



63.6 
88.3 
82.3 
72.3 
91.3 
58.6 
82.5 
82.3 
34.5 
69.4 
65.9 
72.0 
72.2 



NOTES: Models are abbreviated to XSPEC syntax: WA, PO and DISKBB as before, BB-blackbody continuum. Model fitted with a cool 1 or hot 2 
DISKBB component. a External absorption column (10 21 atoms cm -2 ), b PO photon index, c PO normalisation (10~ 3 photon cm" 2 s — 1 keV - 1 
at 1 keV), d inner disc temperature (keV). e DISKBB normalisation as before (((Ri n /km)/(D/10 kpc)) 2 cos G), f fraction of the total flux (0.3-10 keV) 
in the PO component, 9 fraction of the total flux (0.3-10 keV) in the DISKBB component, h blackbody temperature (keV), ' blackbody normalisation, 
(10" s L 39 /D 2 ; where L39 -source luminosity in 10 39 erg s Dio-source distance in 10 kpc). J Fraction of the total flux (0.3-10 keV) in the 
BB component. k Unconstrained at the 90 per cent confidence level. In the three cases in the central portion of the table where no fit is shown, the 
X 2 -minimisation always found the minimum describing the IMBH model. Statistically acceptable fits are highlighted in bold. (NB. In cases where the 
spectral fitting was not performed over the 0.3-10 keV energy range, we created a dummy response matrix in XSPEC (see Sec. 4.5), to determine the 
fraction of the total flux in each spectral component over 0.3-10 keV, consistent with the other datasets). 
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ing to the F-test). However, these were the clearest-cut cases. Of 
the ambiguous spectra, three out of six showed evidence for curva- 
ture. This was marginal in the case of M83 ULX (> 2a level), but 
more significant for NGC 4559 X-l and NGC 5204 X-l (> 3cr), 
with the latter case having an unacceptable power-law fit. Most in- 
terestingly, though, one of the IMBH-fit sources (NGC 1313 X-l) 
shows evidence for curvature (> 3cr improvement), whilst Ho II X- 
1 also shows a significant break (> 4<r improvement). In total 8/13 
sources show at least marginal evidence for curvature at the high 
energy end of their XMM-Newton spectrum. This is strongly sug- 
gestive that the majority of sources are not dominated by a power- 
law continuum at these energies, such as one might expect to see 
if the X-ray emission arises from the optically thin, hot corona as- 
sumed in the IMBH model. 



4.3 Dual thermal models 

In the above analysis we determine that there is (at least marginal) 
evidence of curvature above 2 keV in the majority of our sources, 
arguing that a power-law continuum is not an adequate description 
of this emission, while the hotter accretion disc present in the al- 
ternate model may be more appropriate. The situation below 2 keV 
is less clear cut. Though the alternate model uses a power-law con- 
tinuum to describe the soft component, this is modified to appear 
curved by absorption. As this power-law emission is difficult to un- 
derstand physically we have therefore explored a new variation of 
the alternate model in which the soft component is also intrinsically 
curved {i.e., it does not rely solely on absorption to produce the ob- 
served curvature). In this model we fitted a cool blackbody (BB) 
continuum to the soft component, and a MCD to the hard compo- 
nent. The physical motivation for the soft component comes from 
the suggestion that optically-thick outflowing winds from BHs ac- 
creting at or above the Eddington limit could explain ultrasoft com- 
ponents in ULXs l Kine & Pounds 2003). Such winds would appear 
as BB continua with temperatures of ~0. 1-0.3 keV. One might then 
observe emission from both the accretion disc and the wind given 
a favourable accretion disc geometry and viewing angle (A. King, 
priv. comm.) 5 . 

The results of applying this new model to the data are shown 
in Table [5] It proved to be the most successful empirical descrip- 
tion of the data, providing statistically acceptable fits to 10 out of 
the 13 sources, with 0.15 keV < kT < 0.3 keV and 0.8 keV 
< fcT in < 2.2 keV (excluding NGC 3628 X-l for which the 
temperatures are much higher than shown in the other sources, 
i.e., kT~0.6 keV and fcT in ~4.3 keV). This model provided the 
only acceptable empirical fit to the ULX in NGC 55, and all six am- 
biguous sources are adequately described by this model. Curiously, 
this is also true for both IMBH sources (though it does provide a 
notably worse fit than the IMBH model for M 81 X-9), and for one 
of the two non-standard model sources (NGC 2403 X-l). However, 
this model could not provide statistically acceptable fits to the other 
non-standard source (M 33 X-8) or to the final two ULXs (Ho II X- 
1 & NGC 4395 X-l, which also do not fit well with the power-law 
+ MCD combination). 

5 iMiller et alj l2004al) argue that cool, optically-thick outflows cannot ex- 
plain the soft excess in ULXs as it would be impossible to form powerful 
enough shocks to produce a luminous non-thermal power-law continuum 
component in such systems. If the hard component is not a power-law, but 
instead originates directly from the accretion disc, then the requirement for 
shocks to be present is removed and their argument is circumvented. 




0.5 1 2 5 to 

Energy (keV) 



Figure 2. EPIC pn count rate spectra and Ax residuals for the MEKAL + PO 
model fit to NGC 4395 X- 1 . Also shown are the individual additive model 
components: solid line-PO, dotted line-MEKAL. 

4.4 Why do Ho II X-l and NGC 4395 X-l not conform? 

As shown in the previous sections, statistically acceptable fits for 
Ho II X-l and NGC 4395 X-l cannot be found using the chosen 
two-component models. Interestingly, these are two of the softest 
sources in the sample which means that their X-ray spectra are rel- 
atively dominated by < 1 keV emission, where it is most sensitive 
to absorption characteristics (cf. Fig.0. The analysis of an RGS 
spectrum of Ho II X- 1 bv lGoadetalT l2006) revealed that its X- 
ray emission is subject to absorption by a medium with a sub-solar 
oxygen abundance. The mis-modelling of absorption, in combina- 
tion with the exceptionally high signal to noise EPIC spectrum of 
Ho II X-l, may be responsible for the lack of a good fit to this 
ULX in our analysis. However, even after correcting for the sub- 
solar abundance absorption acting on Ho II X- 1 . iGoad et al J 1 2006 ) 
still did not find acceptable fits using the canonical power-law + 
MCD model. Instead, a statistically acceptable solution to the spec- 
trum was found through using a more physical model, namely the 
DISKPN + COMPTT combination. We explore the utility of a similar 
model in describing our whole dataset in the next section. 

In the case of NGC 4395 X-l, a further inspection of the resid- 
uals to a simple power-law fit does suggest some structure, notably 
including a smooth 'hump' at ~1 keV, which is unlike the fea- 
tureless X-ray spectra of the other ULXs. Therefore we tried a fit 
including a MEKAL component (nominally representing the emis- 
sion spectrum from hot, collisionally-ionised gas) in addition to a 
power-law continuum to model this spectrum. The result is shown 
in Table|8| This model does indeed provide the best (and only statis- 
tically accceptable) fit to NGC 4395 X-l with X 2 /dof = 397.5/366. 
The model fitted to the pn count rate spectrum is shown in Fig.|2| as 
well as the individual additive model components. From this figure, 
one can clearly see how the MEKAL component effectively models 
the emission hump at ~ 1 keV. Using a fixed solar-abundance ab- 
sorber, we find the temperature of the MEKAL component is ~0.75 
keV, while the photon index is quite steep at F ~ 4. 

MEKAL components have been reported in the spectral fits of a 
small number of ULXs. In fact, two other ULXs within this sample 
have been described thus in previous analyses, name ly Ho II X- 
1 bas ed on joint ROSAT - ASCA fits l Miyaii, Lehmann & Hasinaer 
l200ll) . a nd NGC 4559 X-l based on one of two Chandra obser- 
vations <RobertsetaTl l2004). We have investigated whether such 
MEKAL + power-law models can be applied to the current data for 
these sources, and show the results in Table[8| This model is again 
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Table 6. Power-law spectral fits (2-10 keV) 



WA*WA*PO 



X'Vdof 



NGC 55 ULX 
M33 X-8 
NGC 1313 X-l 
NGC 1313 X-2 
NGC 2403 X-l 
Ho II X-l 
M81 X-9 
NGC 3628 X-l 
NGC 4395 X-l 
NGC 4559 X-l 
NGC 4861 ULX 
NGC 5204 X-l 
M83 ULX 



3.58"* 



3+0.06 
-0.07 

2.59±0.04 

1 70+ ' 05 

'• /u -0.07 
Z - J -0.16 

2 63 +4.95 

z - D:, -0.07 

2.61±0.03 
1.78±0.05 

r,+0.06 

-0.08 
3+0.31 
-0.45 
3+0.08 
-0.10 

2.50±0.38 
1.89±0.09 



1.507 
4.48^ 
2.28 _t 



2.67 



+0.21 



8.90Z 
0.44^ 
0.27 4 



3.57±0.25 
f0.41 
-0.39 
Hi -OS 
-0.04 
7 +0.04 
-0.05 
n n 7 +0.05 
U - SZ -0.06 

2.65±0.10 

1.59±0.10 

0.10±0.01 

0.12±0.04 
24 +0.02 

u - z ^-0.03 

09+ 06 

uuy -0.04 

29+ - 04 
u - z -0.03 

0.17±0.05 



459.3/326 
663.1/548 
263.7/259 

60.2/70 
467.3/340 
805.7/824 
360.4/355 
219.5/211 

33.1/28 
171.3/155 

12.3/8 
160.8/133 

45.5/39 



NOTES: Models are abbreviated to XSPEC syntax: PO 



as before. 



PO photon index, 



PO normalisation 



(10 -3 photon cm -2 s _1 keV -1 at 1 keV). Statistically 
acceptable fits are highlighted in bold. 



Table 7. Broken power-law spectral fits (2-10 keV) 



WA*WA*BKNPO 




p b 

-^brcak 


r 2 c 


A B p d 


X 2 /dof 


Ax 2e 


l-P(F-test) / 


NGC 55 ULX 


3.08±0.11 


3 90 +0.20 
JW -0.18 


5 34+0.53 

J "' -0.41 


9 09+O.26 
"-0.24 


341.0/324 


118.3 


>99.9 


M33 X-8 


2.17±0.09 


3 O4+0.26 


o ,9 + 0.19 

-0.15 


g 07+0.56 


516.1/546 


147.0 


>99.9 


NGC 1313 X-l 


, 55+O.IO 
JJJ_0.09 


4 n 3 +1.37 
^■ yJ -0.53 


2 , fi +0.66 
z - ID -0.25 


n 07+O.O5 
u --"-0.04 


250.9/257 


12.9 


99.8 


NGC 1313 X-2 
















NGC 2403 X-l 


2 07 +0.10 
z - u/ -0.11 


4 00+° 17 
^■ uu -0.15 


4 05+0.33 


0.49±0.05 


325.4/338 


141.9 


>99.9 


Ho II X-l 


2.55±0.04 


S ^|+°' 4: - i 
J - J1 -0.56 


3.08ig;» 


250+0.11 
Z - J -0.12 


786.1/822 


19.6 


>99.9 


M81 X-9 
















NGC 3628 X-l 
















NGC 4395 X-l 


9 77 + 1 .31 
z - "-2.58 


2 4 5 +0.93 
z -^ J -0.23 


402+0.77 


<0.09 


30.0/26 


3.1 


72.1 


NGC 4559 X-l 


9 07+0.12 
z - u/ -0.14 


4 77+0-65 
^■"-0.41 


3 14+O.62 


0.20±0.03 


157.7/153 


13.6 


99.8 


NGC 4861 ULX 
















NGC 5204 X-l 


1 65+ ' 10 


4 Q2+ 46 
^■ yz -0.41 


9 qq+0 65 
zyy -0.45 


n 90+0.03 

u - ZJ -0.04 


138.3/131 


22.5 


>99.9 


M83ULX 


<1.79 


2 41+0.31 
z -^ 1 -0.15 


, ,+0.46 
- 31 -0.40 


< 0.08 


37.5/37 


8.0 


97.2 



NOTES: Models are abbreviated to XSPEC syntax: BKNPO-broken power-law model. a Photon index below the 
break energy, 6 break point for the energy (keV), c photon index above the break energy, d BKNPO normalisation 
(10 — 3 photon cm -2 s — 1 keV -1 at 1 keV), e x 2 improvement over a single PO fit, for two extra degrees of free- 
dom, ' statistical probability (per cent) of the fit improvement over a single PO fit. We only show the results where we 
could constrain fits that showed the requisite behaviour, i.e. T2 > Ti and inbreak constrained in the 2-10 keV range.. 
Statistically acceptable fits are highlighted in bold. 



Table 8. MEKAL + power-law spectral fits 



WA*WA*(MEKAL+PO) 




kT b 


A M C 






X 2 /dof 


fv f 


fV 9 


Ho II X-l 

NGC 4395 X-l 
NGC 4559 X-l 


1.42±0.04 
1.73±0.13 
1 26+ 010 

1 - ZD -0.11 


0.66±0.03 
0.75±0.04 
0.25±0.02 


1.32±0.15 

7 90+1.56 
'■"-1.54 


2.67±0.02 

oo+0. 08 
j - SS -0.09 
9 9S-+0.05 
z - ZJ -0.04 


2.90±0.05 

<r on+0.46 

3-3»_ 014 

0.24±0.01 


1463.9/1312 
397.5/366 
612.1/597 


96.7 
84.4 
92.8 


3.3 
15.6 
7.2 



NOTES: Models are abbreviated to XSPEC syntax: WA and PO as before, MEKAL-thermal plasma model. a External absorption column 
(10 21 atoms cm' 2 ), b plasma temperature (keV), c MEKAL normalisation ((10 _19 /(47r[D(l + z)] 2 )) J n c nndV, where D is the 
distance to the source (cm), n c is the electron density (cm -3 ), njr is the hydrogen density (cm -3 )), d PO photon index, e PO normal- 
isation (10 — 5 photon cm -2 s — 1 keV - 1 at 1 keV), f fraction of the total flux (0.3-10 keV) in the PO component, 9 fraction of the 
total flux (0.3-10 keV) in the MEKAL component. Statistically acceptable fits are highlighted in bold. (NB For Ho II X-l we created a 
dummy response matrix over the 0.3-10 keV energy range, to determine the fractional flux in each spectral component (see Sec 4.5)). 
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rejected for Ho II X-l, though it is notable that the quality of the fit 
is very close to that of the best empirical description of the source 
(cool MCD plus power-law). In this case the temperature of the 
MEKAL is again relatively high at ~ 0.66 keV, but the contribu- 
tion of this component is minimal (~ 3 per cent of the 0.3-10 keV 
flux) . This is both far hotter and far fainter than the ~ 0.3 keV 
thermal plasma that comp osed 20 - 30 pe r cent of the 0.5-2 keV 
emission modelled bv lMivaii et alJ 1200 ll) . On the other hand, the 
temperature of the best-fitting MEKAL to NGC 4559 X-l is lower, 
at ~ 0.25 keV, but it does constitute part of a statistically-acceptable 
fit to the da ta. This temperature is similar to the ~ 0. 1 8 keV plasma 
inferred by Rob erts et alJ 12004). We note that this latter result just 
adds to the spectral ambiguity already seen in the case of this ULX. 

Finally, if NGC 4395 X-l really does possess an X-ray line- 
emitting component, what is its physical origin? Several ideas have 
been put fo rward to exp l ain su ch a component in the spectrum 
of an ULX. IMivaii et all {2001) suggest the presence of a young 
supernova remnant coincident with Ho II X-l, a theme expanded 
on bv lFeng & Kaaretl 120051) . Indeed, these authors independently 
confirm the presence of the MEKAL component in the spectrum of 
NGC 4395 X-l, and fi nd two further ULX s with similar character- 
istics. Alternately, as iRoberts et all 120041) find the MEKAL compo- 
nent in NGC 4559 X-l to switch on between two Chandra obser- 
vations separated by ~ 5 months, they suggest that the plasma may 
originate in the collision of a jet or outflow from the ULX with 
a denser medium in the close proximity of the syst em. A second 
alternative is offered by Terashima & Wilson 1 2004), who suggest 
that emission lines in the spectrum of an ULX in M51 may orig- 
inate from a photoionised stellar wind, similar to what is seen in 
some high-mass X-ray binaries in our own Galaxy. A final point of 
interest is that a broad ~ 1 keV feature, similar to that driving the 
MEKAL fit in NGC 4395 X-l, has recentl y been seen in the XMM- 
Newton spectrum of GRS 1915+105 by iMartocchia et alJ 1 2005) 
(although see their paper for caveats). One intriguing possibility 
is that this feature could originate in a disc wind. Clearly ULXs 
possessing MEKAL components in their spectra are an interesting 
subject in their own right, and should be the subject of future atten- 
tion. 



4.5 Physical models 

The success of t he COMPTT mode l in describing th e spectra of 
some ULXs (e.g. iMiller etai]|2003l ISoad et alJl2006l) shows that 
a Comptonised spectrum is in principle another viable alternative. 
The COMPTT model however, is not fully self-consistent e.g., it 
allows spectra to have temperatures higher than expected when 
Compton cooling is taken into account. We therefore attempted to 
fit a m ore physical m odel to the ULX spectra, namely the EQPAIR 
model (CoDoi 1999), which allows thermal and non thermal elec- 
tron distributions. In our modelling, we assume a purely thermal 
electron distribution, the temperature of which is computed self- 
consistently by balancing heating and cooling (the latter of which 
is mainly due to Compton cooling). A key parameter of this model 
is the ratio liJL, where and l s represent the compactness of the 
electrons and the compactness of the seed photon distribution re- 
spectively. This ratio depends on the geometry of the source as well 

6 The contribution of the MEKAL is so small that its parameterisation is 
relatively insensitive to changes in metallicity. For example, setting a low 
abundance (as found by IMivaii et al . 2001) leaves the temperature of the 
MEKAL unchanged at ~ 0.66 keV. 



as the internal mechanisms, and it provides the dominant influence 
on the spectral shape. The source geometry assumed is either spher- 
ical or a disc-corona (slab) geometry. Low energy (UV or X-ray) 
thermal photons from the accretion disc are assumed to be emit- 
ted uniformly inside the source region for the spherical models and 
enter from the base of the corona in the slab geometry. We take 
the seed photons to have an accretion disc spectrum as described 
with the DISKPN model of Gierlihski et al. ( 1999), with the inner 
edge of the disc at QGAl/c 2 . This model assumes a proper general 
relativistic potential and has a characteristic temperature fcTbb- To 
allow for a patchy Compton corona, we add a second DISKPN com- 
ponent to the model, with the temperature coupled to that of the 
EQPAIR seed photons. 

This model ultimately proved the best description of the ULX 
spectra in our sample, with statistically acceptable fits to 11/13 
datasets, all with similar or improved goodness of fits compared to 
the best empirical modelling, for only one extra degree of freedom. 
Only Ho II X-l and NGC 4395 X-l were rejected at > 95 per cent 
confidence (see Table[5), although in the case of Ho II X-l it does 
provide clearly the best fit of all that were attempted to the data, 
and is only marginally rejected (P re j « 97.7 per cent). Again the 
fit to NGC 4395 X-l is poor, for reasons explained in the previous 
section. 

The fits provide a uniformly low measurement of the disc tem- 
perature across the sample, in the range 0.07 < fcT max < 0.29. 
Taken in isolation, this might be interpreted as strong support for 
the presence of IMBHs in ULXs. However, many of the fits show a 
second remakable characteristic, which is that the optical depth of 
the coronae appear to be very high, ranging from r ~ 8 for Ho II 
X-l up to depths well in excess of 30 (NGC 55 ULX, NGC 2403 X- 
l) 7 . Importantly, this result gives a physical explanation for the cur- 
vature noted to be present in the spectra in previous sections. This 
second characteristic appears irreconcilable with the main assump- 
tion behind the IMBH mo del, which is that they operate as simple 
scaled-up BHBs (see also lGoad et alj2 006). as such sources typi- 
cally do not possess very optically-thick coronae. This result must 
therefore strongly challenge the IMBH interpretation for the nine 
ULXs (including Ho II X-l) well-described by this model. How- 
ever, three ULXs - NGC 1313 X-l, M81 X-9 and (possibly) NGC 
3628 X-l - may still possess optically-thin coronae, and as such 
remain viable IMBH candidates. 

As this spectral model provides at least the same goodness of 
fit as the empirical models in all but one case, we used these fits to 
determine the observed X-ray flux. Due to uncertainties in the low 
energy spectrum of some ULXs, the spectral fitting was restricted 
to 0.5-10 keV (and in the case of Ho II X-l the pn data was re- 
stricted to 0.7-10 keV) instead of the 0.3-10 keV range adopted 
for the other sources. Therefore to determine the absorbed X-ray 
flux of each ULX over the same energy range, we created a dummy 
response matrix in XSPEC over the 0.3-10 keV energy range for 
those cases with calibration uncertainties. This dummy response 
temporarily supersedes the response matrix used in the spectral fit- 
ting of those sources, allowing us to examine the behaviour of the 
DISKPN+EQPAIR model over this energy range. The three EPIC 
cameras are consistent to within 20 per cent in 12/13 cases, there- 
fore we derived an absorbed flux (0.3-10 keV) for each of these 
ULXs based on an average of the three measurements. We con- 



7 Indeed, the optical depth for four sources hit the artificial upper limit of 
r = 100 during spectral fitting, therefore we quote values as > 'lower limit' 
in these cases. 
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Table 9. diskpn+eqpair spectral fits 



WA*WA*(DISKPN+EQPAIR) TVjj " fcT max b A D C l h /l s d T P e A E ? X 2 ^of 



NGC 55 ULX 


7 ,1+11.12 

z - lz -o.n 


0.24±0.01 


, -17+0. 50 
1 --"-0.36 


4 54+1-90 
qo '*_2.44 


> 37.2 


0Q 6 +3.94 
u - ?D -0.24 


889.3/827 


M33 X-8 


1 03+ ' 09 
lu;, -0.13 


0.08±° ° 7 


< 4.37 x 10 -10 


o 77+O.48 
J - zz -0.21 




1083l% 859 


1188.0/1112 


NGC 1313 X-l 


2 09+ ' 21 
zuy -0.39 


0.21±H 3 2 


< 0.44 


7 77+2.48 
z - zz _0.46 


7+0.4 

u - z -o.i 


9.99+1I, 2 


658.7/672 


NGC 1313 X-2 


1 44+0-30 


Q26 +0.09 
UZD -0.08 


< 0.27 


, S O + 0.95 

l-83_ 40 


n 2+10-7 
— 7.7 


ft 77+2-80 
u -' z -0.34 


254.1/254 


INLrC A-l 


7 ,-0+0.23 


n 97 +0.05 
°- 27 -0.04 


°- n -o.06 


A sn +7.39 
4 - 50 -1.83 


> OO.O 


n o fi +l-21 
0.36_ 31 


fill 


HoIIX-l 


, n +0.05 
1 - 11 -0.07 


?0+ ' 02 
uzu -0.01 


20 ,+0.27 
Z - UI -0.19 


1.01±0.03 


c+0.6 


26 7+ 4 ' 4 
ZD -'-8.5 


1414.8/1311 


M81 X-9 


1 90+ ' 14 

1;nJ -0.19 


02 3+o.oi 

u - z -0.03 


< 2.08 


o 9s +2.09 
J - z -0.47 


5 +o.3 


lg g+16.2 


873.8/861 


NGC 3628 X-l 


4 28+ ' 36 
^■ zo -0.42 


007+0.02 
uu/ -0.01 


, 6 7+972 
1D '-112 


26 9+ 12 6 
ZD ' y -9.4 


< 7.4 


47 8+ 21 ' 1 
^'■°-11.2 


424.9/403 


NGC 4395 X-l 


1 00+° 12 

luu -0.08 


0.20±0.01 


15+0-08 


, O7+0.34 
J -="-3.17 


> 30.4 


03+0 03 
u - UJ -o.oi 


467.0/365 


NGC 4559 X-l 


1 Q7+0.15 


16+ 001 

uao -0.02 


7 71+3.57 
z - z -0.76 


7 70+O.5I 
Z - ZJ -0.71 


13 n+3.4 


4 12+ 1 ' 49 
^■ lz -1.45 


577.5/596 


NGC 4861 ULX 


1 09+ ' 81 


0.24±0.06 


06+ 22 

u - uo -0.03 


5.30±S;f| 


> 15.8 


06+ ' 55 
u - UD -0.04 


78.7/71 


NGC 5204 X-l 


0.24±0.08 


7Q+0.03 
uzy -0.05 


09+O- 09 
u ' uy -0.05 


4 07+3.27 
^""-1.25 


26 3+ 57 
za;) -5.6 


O25+O.6O 
"■ Z J_0.13 


505.8/526 


M83ULX 


65+ 34 


n 70+0. 08 
u - z:> -0.07 


05+ ' 21 


1 96+ 7 ' 63 


> 15.8 


32+3.98 
U - 5Z _0.19 


193.2/204 



NOTES: Models are abbreviated to XSPEC syntax: WA as before, DISKPN-accretion disc model, EQPAIR-Comptonisation model, 
a External absorption column (10 21 atoms cm -2 ), b maximum temperature in the accretion disc (keV), c DISKPN normalisation 
(10 — 3(M 2 cos(i))/(D 2 /3 4 ); where M-central mass (Mq), D-distance to the source (kpc), i-inclination angle of the disc, /3-colour/effective 
temperature ratio), d ratio between the compactness of the electrons and the compactness of the seed photon distribution, e optical depth, ' EQ- 
PAIR normalisation (corresponding to the disc component) i.e., (f c M 2 cos(i)) / (D 2 [3 4 ), where f c is the covering factor. Statistically acceptable 
fits are highlighted in bold. 



verted these fluxes to an observed X-ray luminosity for each source 
in the 0.3-10 keV band (assuming the appropriate distance) and 
tabulate the results in Table^. 



5 A COMMENT ON THE RELATIONSHIP BETWEEN 
X-RAY LUMINOSITY AND INNER-DISC 
TEMPERATURE IN ULXS. 

Perhaps the most visually striking evidence for ULXs con taining 
IMBHs was presented by Miller. Fabian, & Miller (2004b) (their 
Figures 1 & 2). These authors selected a sample of ULXs with pub- 
lished estimates of Lx > 10 40 erg s _1 , that require a soft excess 
component (at least at the 3er confidence level) in the low energy 
part of a two component X-ray spectrum. They then compared the 
inferred disc temperatures (taken from the soft component) and un- 
absorbed luminosities (0.5-10 keV) of their ULXs to those of a 
number of Galactic BHBs, and found that these ULXs and stellar 
mass BHBs occupy distinct regions of a Lx - kT diagram. More 
specifically, the ULXs are more luminous but have cooler thermal 
disc components than standard stellar-mass BHBs, consistent with 
the ULXs harbouring IMBHs. 

Following on from this work, we have reproduced Lx - kT 
diagrams, based on the inferred unabsorbed 0.5-10 keV X-ray lu- 
minosities for ULXs in this sample (Fig. However, we do this 
for two cases: (i) disc temperatures measured from the IMBH (cool 
disc + hard power-law) fits; and (ii) disc temperatures taken from 
the dual thermal (cool blackbody + warm disc) fits. In both cases 
we only show sources where the spectral fits are statistically ac- 
ceptable, and the fit parameters are well constrained. For this rea- 

8 Although the power-law + MEKAL spectral model statistically provided 
the best fit to NGC 4395 X- 1 , the measured source flux derived from this 
model and the DISKPN+EQPAIR model were in close agreement. For con- 
sistency, we quote the value from the DISKPN+EQPAIR model. However, as 
the ULX lay on a hot column in the pn data we simply use the average of 
the two MOS flux measurements to estimate its luminosity. 



son, we excluded the values from the IMBH fit to M83 ULX, as the 
disc temperature is unconstrained in this case. Additionally, NGC 
3628 X- 1 does not appear on either plot as its unusual disc tempera- 
tures place it outside the range displayed in both cases. In each case 
the luminosity was derived from an average of the unabsorbed flux 
(0.5-10 keV) from each detector. The errors on the disc temperature 
are 90 per cent confidence errors, likewise the luminosity errors are 
based on the 90 per cent confidence errors in the average measured 
flux, as calculated by XSPEC. Finally, we illustrate th e regions of 
parame ter space occupied by the BHBs and ULXs in Mil ler et alJ 
(2004b) by ellipses, with the ellipse representing IMBH-candidate 
ULXs at the upper-left. 

The two Lx - kT diagrams in Fig.Htell very different stories. 
As would be expect ed, our IMBH model results reproduce those of 
lMiileretaljj2004bl) . i.e., the ULXs occupy a very different region 
of parameter space to the stellar-mass BHBs, consistent with larger 
black holes and hence more luminous, cooler accretion discs. How- 
ever, the dual thermal model disc temperatures suggest an alterna- 
tive interpretation. In this case, the ULXs appear to be a direct, high 
luminosity extension of the BHB class, and follow the Lx oc T 4 
trend one would expect from standard accretion discs. Obviously 
one cannot conclude which is the correct interpretation on the ba- 
sis of these plots. However, it does demonstrate very clearly that 
it is the choice of empirical model used to determine the charac- 
teristics of ULX spectra that governs whether one concludes that 
some ULXs contain IMBHs, or whether they are more similar to 
stellar-mass BHBs. 



6 DISCUSSION 

In this work we have analysed the X-ray spectra of a small sample 
of 13 ULXs, that constituted the highest quality (i.e. most photon- 
rich) datasets in the XMM-Newton data holdings as of December 
2004. We demonstrate that the superior collecting area of XMM- 
Newton over its rivals leads to better spectral definition, primarily 
through the rejection of simple single-component spectral models 
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Figure 3. Unabsorbed X-ray luminosity (0.5-10 keV) plotted against ac- 
cretion disc temperatures inferred from the X-ray spectral fits (as per 
iMilleretalf 2 004b). The ellipses represent the regions occupied by ULXs 
(small ellipse) and Galactic BHCs (large ellipse) presented in Figure 2 of 
iMilleretalJ t2004b). The values of the data points displayed in these fig- 
ures are derived from the IMBH spectral fits (top panel) or the dual ther- 
mal model fits (bottom panel), for sources with statistically acceptable and 
well constrained fits for that model. (NOTE: IMBH data points neglect M83 
ULX, which has an unconstrained temperature, and NGC 3628 X- 1 which 
has a lower temperature than the other sources [~80 eV]; dual thermal data 
points also neglect NGC 3628 X-l which has a much higher temperature 
than the other sources [~4 keV]) 



(power-law, MCD) that have adequately described previous spec- 
tra of many of these ULXs. However, we are still photon-limited 
in many of these spectra, to the extent that in six cases we are un- 
able to tell whether a combination of a MCD plus a power-law 
uniquely fits the data with the MCD component at the low or the 
high energy end of the X-ray spectrum. In each of these six ambigu- 



ous spectra, we accumulated no more than ~ 18000 EPIC counts 
(pn and MOS combined). In each spectrum that we resolved this 
simple ambiguity, or excluded both cases as unacceptable fits, we 
had > 20000 EPIC counts (excluding NGC 4395 X-l, as this had 
both an atypical spectrum and an underestimated pn count rate due 
to its location on a hot column). Even then, the introduction of a 
model composed of two thermal components adds an extra layer of 
ambiguity to many of the datasets. This suggests that one requires 
very photon-rich datasets to make progress in descriptive empirical 
XMM-Newton ULX spectroscopy, equivalent to at least 100-ks of 
flare-free data for a 0.2 count s _1 on-axis ULX observed by the 
EPIC cameras, and spectral results from poorer quality data should 
be regarded with caution. 

Issues of spectral ambiguity aside, perhaps our most interest- 
ing result is the probable detection of curvature in the high-energy 
part of the XMM-Newton EPIC spectrum in more than half of the 
sample . Though this has p r eviously been detected in individual 
ULXs iStobbart etai]|2004l; [Foschini et ai]|2004l). or as minority 
populations in ULX samples (Feng & Kaaret 2005), this is the first 
suggestion that the majority of ULXs, across the whole range of 
ULX luminosity, might appear thus. Indeed, we demonstrate that 
the most successful empirical modelling of our sample is achieved 
using a model in which the high-energy portion of the spectrum is 
described by the MCD model, with a soft excess modelled as a cool 
blackbody emitter. 

This obviously adds to the challenges in modelling ULXs as 
IMBHs. Indeed, though 8/13 ULXs present acceptable fits to the 
simple IMBH model of a cool disc plus a hard pow er-law c ontin- 
uum, w e are again faced with the problem noted by Rober tset alJ 
(2005). Namely we detect a very dominant power-law (> 80 per 
cent of the observed 0.3-10 keV flux) that appears too hard for what 
should be a HS (or perhaps VH) source, if we are to believe that 
the inner-disc temperature can provide an estimate of the BH mass 
(1.6 < F < 2.5 for the IMBH candidates in our sample, compared 
to approxi mately 2.1 < T < 4.8 for the HS, and T > 2.5 for the 
VH state - McClintock & Remillarol2003l) . Curvature above 2 keV 
compounds these problems further, as the IMBH model assumes 
that whilst the bigger black hole results in a cooler accretion disc, 
the corona remains similar to that observed in Galactic BHBs, i.e. 
optically-thin and hence modelled with a power-law. In fact, there 
is no reason to suppose that increasing the mass of the BH should 
alter the state of the corona - timing measurements at least appear 
to scale linearly with BH mass for given spectral states, implying 
the properti es of the state s themselves are invariant with mass (cf. 
iDone & cTe rlihski 2005). But the detection of curvature implies 
that either the state of the coronae is changing as the BH becomes 
more massive, or a different physical process is responsible for the 
emission above 2 keV. 

The dual thermal model uses a MCD to model this curvature, 
with inner disc temperatures in the range 1.1-2.2 keV. As Fig [3] 
shows, this provides a natural explanation of ULXs as a simple ex- 
tension of the behaviour of Galactic BHBs, with more luminous 
and slightly hotter discs. In this case, assuming one can exceed 
the Eddington limit, there is no obvious requirement for IMBHs. 
However, even this simple empirical model has limitations; for 
instance, in order to see both the inner regions of the accretion 
disc and an optically-thick outflow simultaneously, one requires 
a very specific geometry and viewing angle. One may therefore 
have to consider further explanations for the origin of the soft ex- 
cess component; we note that one interesting suggestion, that has 
recently been found to work in some AGN X-ray spectra, is that 
the soft excess could arise in relativistically-blurred atomic fea- 
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tures found in the reflecti on spectrum of a photoionised accretion 
disc ICrummv et all20 06). A further concern for this model is that 
similar spectral shapes have not been commonly found in Galactic 
BHBs, with the exception of slight soft excesses (modelled by the 
extension of the power-law component below the MCD) seen in the 
spectra in LMC X-l and LMC X-3 I H aardt et all200lh . One expla- 
nation could be that soft excesses, below the MCD component in 
Galactic BHB spectra, are not seen due to a combination of high 
absorption columns and lack of detector sensitivity at low X-ray 
energies. However, the columns to a significant minority of BHBs 
are not mu ch in excess of the ~ 10 21 cm" 2 columns inferred for 
our ULXs ( McClintock & Remillard 2003), and many have now 
been observed by XMM-Newton and Chandra without such soft 
excesses being reported. Perhaps this is simply telling us that this 
spectral shape is unique to the very high accretion rates required if 
ULXs contain stellar-mass BHs. 



The model that did best in terms of providing acceptable fits 
to the data was also the only physically self-consistent model we 
used in the analysis, namely the DISKPN+EQPAIR spectral model. 
The best fits provided by this model show, without exception, that 
the accretion disc photons seeding the corona are cool. In fact they 
are cool enough, at < 0.3 keV, to be consistent with the accretion 
disc around an IMBH. However, this model also provides a phys- 
ical explanation for the curvature above 2 keV: unlike in conven- 
tional BHBs, it originates in an optically-thick corona. This com- 
bination of a cool disc and optically-thick corona has already been 
observed in Ho II X-l; Goad et al. 1 2006 ) note the cosmetic simi- 
larities between this spectral model and the model of IZhang et alJ 
(2000) of a three-layered atmospheric structure in the accretion 
discs around BHBs. More specifically, their model includes a warm 
layer (fcT~l-1.5 keV, r~10) between the cool optically-thick ac- 
cretion disc (fcT~0.2-0.5 keV) and the hot optically-thin corona 
(fcT~100 keV, r~l) of BHBs, which is responsible for the dom- 
inant component below ~10 keV (see also iNavakshin & Melial 
ll997llMisraetalll998l) . The warm layer appears relatively stable 
and hence unconnected to the hot coron a, which can be highly dy- 
namic and even disappear completely iM isra et alJll99a) . In this 
picture, the cool disc seeds the warm optically thick scattering 
medium, and as such may explain both components seen in our 
modelling. 



This model has been successful in descr ibing the X-ray spec- 
trum o f GRS 1915+105 IZhang et alj|2000h . and so iGoad et ail 
( 2006) speculate that this modelling, combined with an upper limit 
on the mass of the BH in Ho II X-l of ~ 100- Mq derived from 
timing properties, implies that Ho II X-l behaves in an analogous 
manner to GRS 1915+105 in its x-d ass - We speculate, by exten- 
sion, that many of the ULXs in our sample that are well fitted by 
this spectral model could also be GRS 1915+105 analogues, i.e. 
stellar-mass BHs accreting at around and in excess of the Edding- 
ton limit. We note that if individual BHs with m ass up to 80-Mp, 
can fo rm from stellar processes as suggested bv lBelczvnski et all 
(2004), and these larger stellar-mass BHs exist in ULXs, then the 
factors by which the Eddington limit are exceeded need not be large 
for even the brighter ULXs (factors ~ 2 — 3 only). One piece of 
additional, encouraging evidence in this respect is that, similarly 
to Ho II X-l and GRS 1915+105 in its analogous x-class of be- 
haviour, the light curves of the ULXs in this sample show little or 
no intrinsic variability on timescales of minutes to hours, with the 



fractional variability (in excess of counting noise) limited to < 10 
per cent in all cases 9 . 

Recently a second, very plausible explanation for our 
physi cally self-c onsistent spectral modelling has come to light. 
iDone & Kubotal 120051) describe spectral modelling of the Galac- 
tic BHB XTE J1550-564 in its high-luminosity VH state, using a 
model in which the energetics of the inner regions of the accretion 
disc are coupled to a surrounding corona. This results in a cooler 
apparent disc temperature, as the corona drains energy from the 
inner disc, and an optically-thick corona, that are both part of the 
accretion flow. Done & Kubota (2005) note that this may provide 
an explanation for the low disc temperatures observed in ULXs not 
reliant upon the presence of an IMBH. As we also provide evidence 
that the corona itself is indeed optically-thick in such sources, this 
model must constitute a very serious physical alternative to IMBHs 
for the maj ority of ULXs. H owever, even here there are caveats, 
for instance iFeng & Kaaretl 120051) note that optically-thick coro- 
nae should show deep Fe K absorption edges, unless the accreting 
material has a very low metallicity. We do not detect such features 
in our analysis. 

Our empirical spectral fitting has posed new challenges for 
the IMBH model fits to ULX spectra, and our physical mod- 
elling describes scenarios in which the bulk of ULXs could 
be stellar-mass BHs accreting at around the Eddington limit. 
This is perhaps not surprising, as many strands of recent ev- 
idence have pointed away from an IMBH model for most 
ULXs. In fact, excepting the somewhat unique case for M82 
X-l as an IMBH, possibly formed in the dense MGG-11 clus- 
ter or captured as the nucleus of an accreted dwarf galaxy 
(e.g. IS trohmayer & Mushotzkv 2003; Portegies Zwart et alJl2004t 
King & Dehnenl l2005: Muc ciarelli et al]l2006l) . and the cool disc 
detections (shown to be somewhat ambiguous in this paper), ob- 
servational results have tended to argue in the opposite sense. For 
example, the probable breakin g of the Eddington li mit seen in 
some - perhaps most (see Jonker & Nelemans 2004) - Galactic 
BHBs, and especially GRS 1915+105, argues that we cannot ex- 
clude stellar-mass BHs from producing ULXs on this trivial basis 
iMcClintock & Remillardl2003l) . Also, the shape of the universal 
X-ray Luminosit y Funct i on (X LF) for high-mass X-ray binaries, 
derived by Gri mm et alJ 120031) . is somewhat puzzling if IMBHs 
constitute a significant part of the ULX population. In particular, 
why does the XLF appear to cut-off at ~ 2 x 10 40 erg s -1 ? 10 
Surely no such cut-off would be present if ~ 1000-A/q IMBHs 
constitute a large fraction of the ULX population (we are certainly 
unaware of any other source population that cuts off at ~ 0.1 Ed- 
dington rate). We therefore conclude that, on current evidence, it is 
unlikely that accreting IMBHs constitute a large proportion of the 
total ULX population. 



7 CONCLUSIONS 

We have conducted a detailed examination of the X-ray spectral 
shapes in a sample of the highest quality XMM-Newton EPIC ULX 

9 Excepting the N GC 55 ULX, whic h shows prominent dips in its XMM- 
Newton light curve IStobbart et all2004l) . 

10 A near ident i cal cu t-off is found from a large sample of ULXs studied 
by ISwartz et alJ j2004h . Independent support for this cut-of f comes from 
the empirical Lx — star formation rate relationship of Grim m"et all I20O3I) . 
which can o nly have its l inear form above ~ IOMq yr _1 if the cut-off is 
real (see also lGilfanov et all2004l) . 
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datasets available to us. Most notably, more than half of the ULXs 
show at least marginal evidence for curvature in their 2-10 keV 
spectra, which is somewhat unexpected if they are to be interpreted 
as the accreting ~ 1000-Mq IMBHs suggested by modelling the 
soft spectral components as accretion discs. Physical modelling 
shows that this curvature is likely to originate in optically-thick 
coronae, which in turn leads to interpretations of the ULXs in terms 
of high accretion-rate stellar-mass (or slightly larger) BHs operat- 
ing at around the Eddington limit. However, while we conclude that 
it is likely that the general ULX population does not have a large 
contribution from IMBHs, we obviously cannot rule out the possi- 
bility that some ULXs do possess IMBHs. Perhaps the best candi- 
date on the basis of our spectral fitting is M8 1 X-9, which is well fit- 
ted by cool disc plus power-law/optically-thin corona models, and 
does not show explicit curvature in its 2-10 keV spectrum, though 
even this ULX may be fitted using a hot (~ 2.2 keV) accretion disc 
plus soft excess model. Clearly, it is difficult to find unique solu- 
tions for these sources even with high quality XMM-Newton EPIC 
data. Ultimately, we may perhaps have to wait for radial velocity 
measurements from the optical counterpart of an ULX, leading to 
dynamical mass measurements of the compact accretor, before we 
have conclusive evidence whether any individual ULX does har- 
bour an IMBH. 
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